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We study the effects of intrinsic sequence-dependent curvature for a two dimensional semiflexible 
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subjected to any external force, such a system is equivalent to a system with a well-defined intrinsic 
curvature and a proper renormalized persistence length. We find the exact expression for the 
distribution function of the equivalent system. However, we show that such an equivalent system 
does not always exist for the polymer subjected to an external force. We find that under an external 
force, the effect of sequence-disorder depends upon the averaging order, the degree of disorder, and 
the experimental conditions, such as the boundary conditions. Furthermore, a short to moderate 
length biopolymer may be much softer or has a smaller apparent persistent length than what would 
be expected from the "equivalent system" . Moreover, under a strong stretching force and for a long 
biopolymer, the sequence-disorder is immaterial for elasticity. Finally, the effect of sequence-disorder 
may depend upon the quantity considered. 
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I. INTRODUCTION 

It is known that sequence-dependent properties of 
biopolymers play a crucial role in many biological pro- 
cesses. More specifically, sequence-disorder has impor- 
tant influences on DNA packaging, transcription, repli- 
cation, recombination, and repair processes [l], H H 0, 
H, H 0, H, S El EH ■ Owing to progress in experimen- 
tal techniques such as laser or magnetic tweezers and 
atomic force microscopy, it is now possible to manipulate 
and observe single biomolecules directly, and thus make 
a better comparison between theoretical predictions and 
experimental observations. 

In theoretical studies, a semiflexible biopol ym er is of- 
ten modeled as a filamentiJE i, 0, i, El tH El El, 
El, El El E3, El El HI • For instance, the worm- 

like chain (WLC) model, which views the biopolymer 
as an inextensible chain with a uniform bending rigid- 
ity but with a negligible cross section, has been used 
successfully to describe the entropic elasticityof a long 
double-stranded DNA (dsDNA) H, El, El, El ■ How- 
ever, the traditional elastic models are usually uniform 
and ignore the role of sequence-disorder. Under what 
conditions such a simplification is valid is therefore an 
intriguing question. Based on the elastic models, two ef- 
fects of sequence-disorder need to be considered. First, 
structural inhomogeneity yields variations of the bending 
rigidity along the filament, and results in an s-dependent 
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persistence length l p (s) @, [1], where s is the arc length. 
It has been reported that for a long biopolymer with 
short-range correlation (SRC) in l p (s) and free of ex- 
ternal force, this effect can be accounted by a replace- 
ment of the l p (s) by an appropriate average HQ. How- 
ever, for a short biopolymer, inhomogeneity in l p (s) tends 
to make physical observables divergent [1, Q. Second, 
the local structure can be characterized by the intrinsic 
sequence-dependent curvaturesfi.e.jthe static curvature 
or the frozen-in curvature) @, H II [1, S 0, 1, 1 El EH , 
and this is also the focus of the present work. For a 
long biopolymer with zero mean curvature, again it has 
been demonstrated that the effect of intrinsic sequence- 
dependent curvatures can also be reduced into a simple 
correction of the uniform persistence length, either free 
of external force [1, d, [H [2J| or under moderate external 
force [H EH, Hl| . However, it is argued that sequence dis- 
order is immaterial for the elasticity of a long DNA under 
strong stretching force [H, EH- In this paper, we prove 
it exactly for a two dimensional biopolymer even with a 
nonvanishing mean intrinsic curvature. Moreover, it is 
well known that the short or intermediate-length DNAs 
play a more important role than the long DNAs in bio- 
logical processes, from DNA packaging, to transcription, 
gene regulation and viral packaging [H [H [Hj]. As a 
consequence, the effect of intrinsic sequence-dependent 
curvatures for short or intermediate-length biopolymers 
requires more attention. When the biopolymer is free 
of external force and with SRC in intrinsic sequence- 
dependent curvatures, it has been shown exactly that 
such a three-dimensional (3D) system is equivalent to 
a system (we will refer it as the "equivalent system" 
henceforth) with a well-defined (i.e., without random- 
ness) intrinsic mean curvature and a corrected persis- 
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tence length [1, irrespective of its length. In this 
work, we show that the same conclusion is also valid 
in the two-dimensional (2D) case and present the gen- 
eral solution of the distribution function of the equiva- 
lent 2D system. On the other hand, the effects of intrin- 
sic sequence-dependent curvatures in a short biopolymcr 
under external force are not yet known. In this paper, 
we demonstrate that under external force, the effect of 
sequence-disorder for a short biopolymer is dependent 
on the average order and the experimental conditions. 
Moreover, we find that the results arc also dependent on 
the boundary conditions (BC), and the short biopolymcr 
looks softer than what is expected from the "equivalent 
system" . Because theoretically a 2D system is relatively 
easier to study and experiments on the conformations of 
biopolymers are often conducted in a 2D environment 
@, Hil i this work will focus on the 2D system. 

The paper is organized as follows. In section II we 
set up our model. Section III presents the exact proof 
of the existence of the "equivalent system" and the gen- 
eral distribution function of the "equivalent system" for 
a force-free biopolymer. In section IV we focus on the 
conformation and elasticity of the biopolymer under con- 
stant external force. Follows a section discussing the ef- 
fects of disorder in a segment dependent curvature in the 
constant extension ensemble. Finally, we end the paper 
with conclusions and discussions. 



WLC model [1J, [13|, |14|, [15( . Note that with free bound- 
ary condition, a negative value of f x only extends the 
polymer in the negative direction rather than the positive 
direction. So it does for a long polymer since the bound- 
ary condition becomes unimportant. Therefore, the sign 
of the force is meaningless in these cases. However, it is 
not the case for a short polymer with a fixed initial angle. 

If both l p and c(s) are well-defined functions of s, a 
macroscopic quantity is defined as the average with 
Boltzmann weights over all possible conformations, 



(4) 



where Zk = / T>[4>]e~~ £ . Function B[{</>(s)}] represents 
different physical situations. For instance, if 2?[{0(s)}] = 
t(si) • t(s2), we find the orientational correlation func- 
tion (OCF); if B[{(f>(s)}] = \r L - r | 2 , we obtain the 
mean end-to-end distance, where vl — r(L) and r = 
r(0); if B[{4>(s)}} = 8(TL - tftds), we get the dis- 
tribution function of end-to-end vector; if £?[{</>(s)}] = 
8{yl — r )S[t(L) — t(0)], we find the looping probability. 
B[{(f>(s)}] is independent of k and c(s) but can be a very 
complex function of <j)(s) and 4>(s). 

Eq. ([I]) uses <fi as the variable of integration. However, 
the variable of integration can be replaced by </>(s), i.e., 
we have the following identity (see Appendix 1 for proof) 



II. THE MODEL 

The configuration of a filament with negligible cross 
section can be described by the tangent vector, t(s), to 
its contour line, where s measures the location along the 
filament. In two dimensions, t(s) = {cos0(s),sin0(s)}, 
where the azimuthal angle <f> is the angle between the 
x-axis and t. The locus of the filament can be found by 



r(s) = {x(s),y(s)} = / t(u)du. 

Jo 



(1) 



The reduced energy of the filament with intrinsic curva- 
ture but free of external force can be written as: 



£ [{<Ks)}] = 



E 
k^T 



c(s)] 2 ds, (2) 



where E is the energy, <j> = d<p/ds, T is the temperature, 
kg is the Boltzmann constant, L is the total arc length 
of the filament and is a constant in the model so that 
the filament is inextensible, k = l p /2 with l p the 2D bare 
persistent length, c(s) is the intrinsic sequence-dependent 
curvature. Under a uniaxial applied force f x (along x- 
axis), the reduced energy of the filament becomes 



£ = £ [{(j)(s)}] - f I ens < >(/.-•. 

'0 



(3) 



where the reduced force is defined by / = fx/ksT. When 
c(s) = and l p is a constant, it returns to the well known 



fV[<i>( S )}B[{cl>(s)}}e- £ fD[<j>( S )]B[{<j>(s)}]e- 



fi>[<K*)]* 



(5) 



For a biopolymer without correlation on c(s), or with 
SRC on c(s) but in the coarse-grained model, the distri- 
bution of c(s), W({c(s)}), can be written as a Gaussian 
distribution with mean c, and root-mean squared devia- 
tion l/Ja; 



W({c(s)}) = exp 



-[c(s) -c] 2 eLs 



(6) 



In this case, we need to average over c for all biopolymers 
in the system. Note that averaging can be done in two 
different orders, cither 

B = (B[{ct>(s)}}) = ±J V[c( S )}W({c( s )})B^, (7) 

where Z a = J V[c(s)]W ({c(s)}) , or 



B' = (B[{<t>(s)}])' 



1 

Zk 
1 

Zk 



1 



V[cj>]B[{4>} 



V[c]W{{c})B[{<j>}]e 



1 



V[c]W({c})e- 



.(8) 



Physically, Eq. ([7]) corresponds to performing a confor- 
mational or thermal average over an individual sample 
first, and then a disorder average over all samples in 
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the system, so we call this a thermal-first system. In 
contrast, in Eq. ([5]) a disorder average over an instan- 
taneous conformation of all samples is first performed, 
followed by a conformational average over all possible 
conformations, so we refer to the corresponding system 
as a disorder-first system. To find a macroscopic quan- 
tity in the disorder-first system is equivalent to averaging 
the sequence disorder over all samples first to create an 
"equivalent" system which is then thermally averaged. 
So it is always possible to construct an "equivalent sys- 
tem" in that case. However, in experiments and compu- 
tational simulations thermal-first averages are performed. 
Therefore, this work will also focus on the thermal-first 
system and show that in many cases the two systems are 
not equivalent. 



III. DISTRIBUTION FUNCTION FOR THE 
FORCE-FREE SYSTEM 

We first present a brief description of the force-free 
system. Using the identities Eqs. §5§ and 



J V[c]W({c})e- £o 



7° 7 

fVf" -Ho 

7° 
H 



(9) 



and exchanging the order of integration, Eq. ([7|) becomes 

JV[c(s)}e- £ W({c(s)}) 



k J 



where 



(10) 



(11) 



Ho = - k[4>(s) - c] 2 ds, 
z Jo 

Zl = J VW~ £ \ Z° = J V[4>]e- H \ (12) 



and the effective persistent length, 

if =2n = 2ka/(k + a). 



(13) 



Note that Eq. (jTUJ) is valid for any L and even if c, k and 
a are s-dependent. Comparing Eqs. ((4]) and (fT0|) . we 
reach the conclusion that a system with SRC in c(s) is 
equivalent to a system with a well-defined mean intrinsic 
curvature c and a renormalized persistence length if. 

The same conclusion has been achieved for 3D biopoly- 
mers following similar arguments @, except that in the 
2D case we have to derive Eq. ([5]) first due to the con- 
vention of using <p a s integral variable. This conclusion 
also means that the thermal-first system is exactly the 
same as the disorder-first system in the force-free case. 

From the standard connection between the path inte- 
gral and the Schrodinger equation, we can find that the 
partition function Z^((f>(s), s; </>(sq), s o) (= / T>[(f)}e~ n °) 



for the system with effective energy H q sa tisfies the fol- 
lowing partial differential equation [2l|, [23, [3l[ 



1 d 2 _ _d_ 



(14) 



Fixing 4>(s) at s = sq, the boundary condition (BC) be- 
comes 



s ; ip, so) = S[(f>-(p], 



(15) 



where ip = 4>(so). In an experiment, one usually takes 
0(0) = when / = 0. 

It is straightforward to show that the normalized func- 
tion (i.e., the distribution function with s > sq) 



P((j),s;cp,So) 



2ttA{s,s ) 



_ e -l<t>-P-f* c(s)ds] 2 /2A(s,s ) 



(16) 



satisfies Eqs. (IT4")) and (fl"5j). where A(s, s ) = f so ds/n(s). 
Eq. (fl6|) can also be derived directly by using a standard 
path integral technique [13, [H| (also see Appendix 2). 



IV. CONFORMATION AND ELASTICITY OF 
THE SYSTEM UNDER EXTERNAL FORCE 

A. On the Disorder-first System 

When / 7^ 0, for the disorder-first system, the deriva- 
tion leading to Eq. (fTU|) can be generalized easily to 
obtain an equivalent system with the effective energy 



H = Ho - f 



is, 



(17) 



no matter what the force may be since the force term in 
£ is independent of c(s) . The equivalent system has been 
well studied 0, HU, [33| . 

In the three dimensional case, we can also reach a sim- 
ilar conclusion by a direct generalization of the proof in 
Ref. Q to find an equivalent system, and an alternative 
proof for the three dimensional filament under weak force 
can be found in Ref. [Ioj |. 



B. General Expressions for the Thermal- first 
System Under Weak Force 

Note that, mathematically, the conclusion of the exis- 
tence of an equivalent system in the force- free case results 
from the fact that both Z® and Z a are Gaussian inte- 
grals so are independent of c(s) or cj). However, such an 
argument fails for the thermal-first system with / 7^ 0, 
because Zk is no longer a Gaussian integral and is depen- 
dent on c, so exchanging the order of integration does not 
simplify the expression. In other words, there is no sim- 
ple way to remove the randomness in c(s) so there is in 
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general no "equivalent system" even under weak force, as 
we will show exactly below. In this case, to first order in 
/, as shown in Appendix 3, 



B = 

with 
B x = 
and 

B 2 = 

where 



^- J V{cj)}Bim}e F [/ V[c]±W{{c})e 



-So 



B\ — B 2 



(18) 



(B) K + f / ds'(B[{<i>(s)}]co S( t>(s')) K , (19) 



/ / d s 'e-^ s ',o)/2 {B [{4>(s)}} costs') + <t>o]) K , 
Jo 

(20) 



A'(s,s ) 



ds/k'(s), k'(s) = ^±|), (21) 
a + 2k 



7 ( S ) = [ S j(s)ds, 7( S ) = M±^£, (22) 
Jo a + k 

and 0o = 0(0) is the initial azimuthal angle. Note that 
j(s) is in general dependent on 0(s) making these ex- 
pressions very complex. 

We should remind that when f ^ 0, (f> is not neces- 
sarily zero but is dependent on the experimental condi- 
tions. In experiments, 4>q may be fixed. In this case, the 
boundary condition at s = is given by Eq. (|15l) . How- 
ever, experiments on stretching biopolymers usually in- 
volve attaching the two ends of the biopolymer to beads, 
and it does not seem to be easy to completely prohibit 
the rotation of the beads. As a consequence, it may be 
difficult to fix (j) . In the extreme case, the polymer can 
rotate freely around the origin. This can be realized by a 
magnetic tweezer [34j . In the more general case, (f>o may 
have a distribution and therefore finally we need to aver- 
age over 4> . It has been reported that different boundary 
conditions have considerable effects on the mechanical re- 
sponse of a homopolymer [l9, 34]. In this work, we come 
to the same conclusion for a heteropolymer. 

On the other hand, for the disorder-first system, we 
find 

B' = ± j V\<t>]B\{<t,}]e- H «Bi- B> 2 , 



B' 



B' = 



f(B[{<Ks)}]) K / d* <«»(#*))), 



f(B[{<P(s)}]) K B' 2 , 
L dse- A ^°V 2 cos \ 



(23) 



c(s)ds 



where Z H = J V[<j>\e- n . 

It is clear that in general B 2 ^ B' 2 since 7 is dependent 
on (p. It in turn leads to in general fi / B'. In other 
words, in general it is impossible to find an equivalent 
system for the thermal-first system under external force. 



Elasticity of the Thermal-first System Under 
Weak Force 



To figure out how serious the effect of the disorder in 
c(s) or how large the discrepancy between B and B', we 
examine the most interesting and also the simplest case 
with constant fc, a and c = 0. It corresponds to the 
WLC model and is often used to describe the entropic 
elasticity of biopolymers, such as dsDNA and proteins. 
Experiments in 3D dsDNA found that k w 78nm and 
k w 45nm[H EI EE Ez|. It follows that k w 1.7k for 
DNA. In this case, 7(s) = k[4>{s) — <fio]/a, and from Eq. 
(f2"3")l . we can obtain 



B' 



iKf (l - c- L / 2K ) cos o (B[{4>(8)}]) K . (24) 



The extension X in the thermal-first system can be 
found as (see Appendix 4) 



X 



(x) - (x) /=0 
Ana 



L _ k{k + 3a) k(k + a) c _ L/k 



k + a k — a 

k 2 cos(20o)e~ 2i / K 7 



L/2k 



X 



a 



3(2fc + a)(3fc + a) 



1 



k- 
6k 2 (1 

+ka (9e 2L/K - 16e 3i/2K + 12e i/fc - 5 

2 ^ e 2L/ K _ 8e 3L/2« + 6e L/fc_ 1 



X, (25) 



On the other hand, in the disorder- first system, the 
extension X' is: 



X' 



2k/ 



" W/=0 

L + k - k (<r L/2K - 2 

2V« _ 6e - L / K + 8e" L / 2 



« 2 / 



cos(20 o ) 
(26) 



From Eqs. (|2"5)) and (|2^|) . we can show that X — > A' 
when a — > 00, as it should be since in this case the system 
is free of randomness. Moreover, for a long polymer, we 
obtain X w X' w 2kL/, so that the averaging order 
is irrelevant for a long polymer, and agrees with known 
results [IE [n|, El]. 

However, the discrepancy between the two systems is 
serious up to moderate length at finite a. From Eqs. 
(1231) and (|2"6"|) , we find that the extensions consist of two 
terms. The first term is independent of BC or </>o, but 
the second term is dependent on cos2</>o so is dependent 
on the BC. As a consequence, different boundary condi- 
tions have strong effects on the elasticity up to moderate 
length. Without loss in generality, we consider two ex- 
treme cases. The first case is to fix </>o = 0, which gives 
the most stringent BC effects. The second is to let 0o 
free so (cos 2</>o) = 0, and only the BC-independent term 
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effect is getting stronger with increasing randomness, i.e. 
with decreasing a. Especially, at k « 1.7k which corre- 
sponds to the experimental value of DNA, we find that X 
can be only about half of X' for L ~ k, and 20% smaller 
than X' up to L ss 10k or about 1500 bp. That means 
DNA is more likely to be in coil state or appears softer 
and has a smaller apparent persistent length than that in 
the "equivalent" system up to a rather long length. Fi- 
nally, these results are independent of the external force, 
so that a small force may produce a considerable effect. 
This fact may be important in a stretching experiment 
for a short polymer since it suggests that the interaction 
between experimental apparatus and polymer may affect 
the results seriously. 



FIG. 1: The ratio r e = X/X' vs L for the fixed BC ( (j>o = 0). 
From top to bottom, the parameters are: k = 1.2k, q = 
6.0K(dash); k = 1.4k, a = 3.5k (dot); k = 1.6k, q = 2.67k 
(dash dot); k = 1.8k, a = 2.25k (solid); k = 2.0k, a = 2.0k 
(short dash). The unit of L is in k. 




FIG. 2: The ratio r e = X/X' vs L for the free BC ((cos 20 o ) = 
0). From top to bottom, the parameters are: k = 1.2k, a — 
6.0k (dash); k = 1.4k, a = 3.5k (dot); k = 1.6k, a = 2.67k 
(dash dot); k = 1.8k, a = 2.25k (solid); k = 2.0k, a = 2.0k 
(short dash). The unit of L is in k. 



remains. For a better comparison, we define a ratio of 
the extensions, r e — X/X' . 

Figures Q] and [2] present some typical results for r e . 
From these two figures, we can see that r e increases 
monotonously with increasing L, and in general the dis- 
order in c(s) makes the biopolymer in the thermal- first 
system softer, i.e. r e < 1, than the "equivalent" system, 
and the effects may be still rather serious up to the length 
of about 20 k. Furthermore, we also find that the results 
are sensitive to BC and randomness of c(s). At first, the 
effect is much more serious (maybe about twice) in the 
free BC (Fig. © than in fixing BC (Fig. [J). Second, the 



D. Orientational Correlation Function (OCF) and 
End-to-end Distance for the Thermal-first System 
Under Weak Force 

Let B[{4>}\ = t(s) • t(s') = cos[<f>(s) - (f>(s% from Eqs. 
(fT5)) - (f2"2")) , we obtain the orientational correlation func- 
tion (s > s') for the thermal-first system and under weak 
external force (see Appendix 5) 



(t(a).t(a')) 



-(s— s ) /2k j_ 



KCOS (f) f 



3(3fc + a)(k - a) 

S = -6k(k - a -) e - 2 l( a + 2k ) s + (a+3k)s']/2ka 

+8a(2fc + a)e- s / 2K - 6fc(3fc + a ) e -[^+(k- a )s']/2k a 



+(k-a)(3k + a) 



2 e -[(a+k)L+(a+3k)(s-s')]/2ka 



-3e" [L+3(s ~ s ' )1/2K + 6e~ s ' /2K - 3c~ (l+s '- s)/2k 
+2e~ (4s-3s ' )/2K + 3 c -[("+ fc ) i +( Q - fc )( s - s ')]/2feaj (27) 

In contrast, in the disorder- first system, 

(t(s) ■ t(s'))' w e -(^ s ')/ 2K + ^k/ costs', 

5' — 6e~ s '/ 2fi — 8e -s / 2K + Qe~ ( - L+S ~ S '^ 2K + 2e^ 3s ' _4s -'/ 2K 

_3 e -[i+3(a-a')]/2« _ 3e -(L+s'- s )/2^ ^8) 

The end-to-end distance can be found as 



with 



AkL 



-L/2k 



4k 2 cos (j>of 



9(fc - a) 2 {k + a)(2k + a){3k + a) 



ry, (29) 



y = l8(k-a) 2 (k + a){2k + a){3k + a) 2 L 

-9k(k - a) 2 (3k + a) 2 (4fc 2 + lla 2 + 22ka) 
-16a(2jfc + a) 2 [3(k - a)(k + a) (3k + a)L 
+Aka(5k + a)(3k - a)}e~ L/2K 
+18ka(k + a) 3 (k - a )2 e -(Vfc+2/a)i 
+18fc(fc + a) 3 (2k + a)(3k + a) 2 e~ L/k 
-ka(k - a) 2 (2k + a)(3k + a) 2 &- 2L ^. 
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On the other hand 



R' 



4kX 



l_^(l_ e -V2» 



jk 2 cos </) /y 

(30) 



with y = 18L - 99k + 16(3L + 4/t)e" 



-L/2k 



+36Ke~ L / K - Ke- 2i/K 



When 



oo, i.e. a system without randomness, we 



obtain i? = i? as it should be. Moreover, L — ► oo also 
leads to i? 2 — > i?' 2 , again supporting the conclusion that 
the order in averaging is irrelevant for a long polymer 

nam eh. 
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FIG. 3: The ratio r fl = i? 2 /-R' 2 vs L - f = 0.5/« and <^ 
( ) (the fixed BC) for all curves. From bottom to top, the 
parameters are: k — 1.2k, a — 6.0k (dash); k = 1.4k, a — 
3.5k (dot); k = 1.6k, a = 2.67k (dash dot); k = 1.8k, a — 
2.25k (solid); k = 2.0k, a = 2.0k (short dash). The unit of L 



Moreover, from Eqs. |gHJ) and |[gDj). we find that the 
end-to-end distance can also be divided into two terms. 
The first term is the force free term and is independent 
of BC, but the second term is dependent on both BC and 
force. Especially, for free BC ((cos2</>o) = 0), the force 
has not effect at all, this is quite different from the force- 
extension relation. Furthermore, with finite (cos20o), R 2 
is also smaller than R' 2 , similar to the force-extension 
relation. However, the ratio tr = R 2 /R' 2 is no longer a 
monotonic function of L, but has a minimum at L w 6k. 
Moreover, the discrepancy between R and R' increases 
with increasing / (see Eqs. (|29l) and (|30p h Note that 
our results are valid only at weak force, and for DNA 
it means that we require ksTf < f c = ksT/n. Taking 
the generally accepted value k sa 50 nm, we have f c « 
0.08pN. Fig. © shows some typical results at / = 0.5/ k, 
which corresponds to an external force of about 0.04 pN. 
From Fig. ([3]), we can see that the discrepancy between 
R 2 and R' 2 is much smaller than that for extensions. 
This is because the force free term dominates the value of 
rft. The fact that the disorder in c(s) has quite different 



effects on the extension and the end-to-end distance may 
be important in experiments. 



E. Elasticity of the Long Biopolymer Under Large 
Force 

Again we assume constant fc, a in this part, but allow 
a nonvanishing c(s). In the large force limit the filament 
is nearly straight, thus t(s) is nearly pointing along the 
direction of the force. That means <p(s) ~ and the 
reduced energy becomes 

£^J^[j>-c{ s )] 2 d s + f^ l -4> 2 ds, (31) 

where we have dropped a constant term —fL. For a very 
long filament, we can use periodic boundary conditions 
with negligible error, so take q n — lixn/L with integer n, 
and expand 4>, c(s) and c(s) as Fourier series 



n=l 



a n sm(q n s), </>(s) 



71=1 



c(s) 



^c„cos(g„s), c(s) = 

n=0 n=0 



a n q n cos(q n s), 

(32) 

c n cos{q n s). (33) 



Note that to use a sine series for cf>(s) is reasonable since 
</>o = <f>(L) = 0. But in general c(L) ^ so we cannot use 
sine series for it. We can also use the full Fourier series for 
c(s) and c(s), but it is straightforward to show that the 
sine part in the full Fourier series make no contribution at 
all so we disregard it. Using the orthogonality property of 
Fourier modes, we can reexpress the energy and extension 
as 



£ 



fcCQ-ZL 



kc^L 



kL x ^ 
T ^ 

n=l 



4 <H 

71 = 1 



{a n q n 



\ 2 f ^ \ ^ 2 

71=1 



{kq 2 n + f)d 2 n 



kf 



Wn + f 




(34) 
(35) 



2 kql + / 



where 



L 



1 



^ OO 



71=1 

£ (fegn) 2 

4^(fcg 2 +/) 2< 



(36) 

(37) 
(38) 
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From Eqs. @ and ([34 )) -((3"8 ]) . we see that in the thermal- 
first system (x) — (xi) + (x 2 ). X\ is independent of c(s), 
and we recover the well known result for the WLC model 



{XI 



1 



1 



(39) 



In contrast, x 2 is independent of 4>(s) and is determined 
by c(s). Rewriting VF({c(s)}) as 



W({c(s)}) 
we obtain 

fa) = 



c -^[(co-e ) 2 + ^2:~ = i(^-5.) 2 
(kq n ) 2 



(40) 



E 



4 ^ (kq n + f) 
L 



<4> 



E 



(fcg») 2 f-2 , 2 

2 K + 



1 (fcg 2 . + ff 



4e 



(fcg„) 2 



1 (fc? 2 + /) 5 



(41) 



On the other hand, replacing k and c„ by k and c n re- 
spectively in Eqs. (|3"4")) - (l3"51) , we find that in the disorder- 
first system the extension becomes 



(x) _ (xi)' + x 2 



L 



L 



= 1 - 



1 



2 x n 



off 



X2_ 

L ' 



(42) 



From Eqs. (|36| - (j42j) . we find that replacing k by n one 
goes from the thermal-first system to the disorder-first 
system. It is also interesting to note that a or the width 
of the distribution of c(s) plays no role in the extension 
in both systems. Furthermore, we find that when c is a 
constant, vanishing or nonvanishing, (x 2 ) ~ since all 
c n = if n > 1. On the hand hand, considering the 
special case with c = ere _As , in the thermal-first system 
we have 



L 



I 00 

iE ( ^ 



(kq n f 



I) 2 " 
{kq n f 



Acrfc 2 



4tt 



q 2 dq 



Acrfc 2 



+ /)2(g2 + A 2) 



(43) 



A long biopolymer in general has a vanishing or small 
mean curvature, or in other words in general a and A are 
small, so that x 2 is negligible in either thermal-first or 
disorder-first systems. Since the corrections of c are neg- 
ligible in both constant and fast decay cases, we can con- 
clude safely that x 2 is always negligible. In the thermal- 
first system, it means that sequence-disorder has no ef- 
fect, which agrees with Marko and Siggia's argument that 



under a strong stretching force, disorder in sequence is 
immaterial for elasticity jiol [l5| . Note that the effect of 
the sequence-disorder has been absorbed into n in the 
disorder-first system, the result in this part provides an- 
other evidence of the non-equivalence of the two systems. 



V. ON THE CONSTANT-EXTENSION 
ENSEMBLE 

Up to now, our discussions are based on a constant 
external force, or in the constant-force ensemble. How- 
ever, the experiments may be performed with fixed ends, 
or in the constant-extension ensemble. For a long poly- 
mer, it is believed that these two ensembles should yield 
the same mechanical properties. However, it has been 
reported that the two ensembles are not always equiva- 
lent for a short polymer. It is therefore interesting to ask 
whether sequence-disorder has the same effects in the two 
ensembles. In the constant-extension ensemble, without 
sequence disorder, the partition function is 



-So 



V[cl>}8(r L - b)e 



= J V[4>}8 ^ tds - b^j 



-So 



(44) 



where b is the constant end-to-end vector. With sequence 
disorder, the disorder-first average for the ensemble be- 
comes 



B' = 



1 

~Z ( . 
1 



v[<p] 



1 



-2a 

V[4>\5(y L - b)S[{<^} 



V[c}W{{c})S{v L -h)e- So B[{4>}} 
V[c]W({c})c- £o 



1 

2 n 



V[4>]6(r L - b)£[{0}]e- 



-Ha 



(45) 



where we have used Eq. §§§ again. Now considering the 
special case -B [{</>}] = 1, which results in B 1 = 1, and 



7° 



1 



Z e ZQ /2?[#( ri -b)e-«o' 
leads us to the result, 



B' = 



JV[mr L -h)Bm} 
fT>[(t>}5{r L -b)e-n° 



-"0 



(46) 



(47) 



Therefore, we reach the conclusion that in the constant- 
extension ensemble for the disorder-first system we can 
still find an "equivalent system" with the effective energy 
Ho. 

On the other hand, the thermal-first average for the 
constant-extension ensemble can be written 



B 



1 

2,. ( 



V[c]W{{c})B e , 
B e = i- y V[<I>]B[{<I>}]6(t l - b)e 



-So 



(48) 
(49) 
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From Eq. (|44]) , we find that due to the existence of the 
8(yl— b) term, Z e is in general dependent on c(s), so that 
an exchange in the order of integration cannot simplify 
the expression for B. In other words, the existence of 
an "equivalent system" for the thermal- first system in 
the constant-extension ensemble is still an open question. 
From our experience in the constant-force ensemble, such 
an "equivalent system" does not exist. 

Our discussions in this section can be directly gener- 
alized to the three dimensional system, so it completes 
and reassess the results of our previous work [8J . 



VI. CONCLUSIONS AND DISCUSSIONS 

In summary, we present a rigorous proof that when 
free of external force, a 2D semiflexible biopolymer with- 
out correlation or with SRC in intrinsic curvatures c(s) is 
equivalent to a system with a well-defined intrinsic cur- 
vature and a renormalized persistence length. We ob- 
tain exact expressions for the distribution function of 
the equivalent system. These conclusions can simplify 
theoretical studies of semiflexible biopolymers, since the 
disorder in c(s) is completely removed in the equivalent 
system. For the system under external force, we find that 
the effect of sequence-disorder is dependent on the order 
in which the averaging is done or the experimental con- 
ditions. In the disorder-first system, it is always possible 
to find an "equivalent system" , no matter the external 
force, the length of the polymer, the statistical ensemble 
or the dimension of the system. However, in the thermal- 
first system, there is in general no "equivalent system" 
for a biopolymer up to moderate length. Physically, this 
is because in the thermal-first system the intrinsic cur- 
vatures favor defects such as kinks, buckles and loops. 
To straighten these defects costs extra energy and there- 
fore requires a larger force. In contrast, the disorder- 
first system erases these extra defects before the appli- 
cation of the force. We find the closed-form expression 
for the force-extension relationship for the elasticity of a 
long biopolymer under a strong stretching force. In the 
thermal-first system, we show exactly that in this case se- 
quence disorder is immaterial even if the biopolymer has 
a nonvanishing mean intrinsic curvature. Moreover, we 
find that in the thermal-first system, the results are also 
dependent on the boundary conditions, and the sequence- 
dependent effects are much more serious in the case of 
free BCs than for fixed BCs. Meanwhile, the results are 
dependent on the degree of randomness and the larger 
the randomness, the more serious the effect. Our results 
suggest that the short biopolymer may be much softer so 
has a smaller apparent persistent length than what the 
"equivalent system" in a mechanical experiment would 
predict. This fact implies that in experiments the in- 
teraction between experimental apparatus and polymer, 
though may be weak, may affect the results seriously for 
a short polymer. Furthermore, our results suggest that 
the effects of sequence-disorder is dependent upon the 



quantity measured and how it is measured. We should 
note that due to the existence of an "equivalent system" , 
the disorder- first system is much simpler in theoretical 
studies. However, it is difficult to realize the disorder- 
first system in experiment. 

On the other hand, we considered weak forces and large 
forces but not intermediate size forces in this work, but 
we can expect that in this case there will also not be an 
"equivalent system" even for a long polymer since this 
is the case for the system under a large force. How the 
elasticity of a system goes from an "equivalent system" 
when free of force to a disorder free system under a large 
force would be an important question to address. We 
also do not consider the system with LRC in c(s), which 
deserves further investigation. 

It should also be noted that this work focussed on 2D 
systems. However, whether the results apply to 3D sys- 
tems is an intriguing question. Mathematically we should 
reach similar conclusions since in the thermal-first system 
exchanging the order of integration does not simplify the 
problem. However, physically there exist some funda- 
mental distinctions between the 3D case and the 2D case. 
At first, the much stronger fluctuations in the 3D system 
may dominate the intrinsic disorder so the effect may be 
suppressed. This would explain why the disorder in c(s) 
has a much smaller effect on the end-to-end distance than 
on the extension. Moreover, the geometry of a polymer 
with natural curvature is also very different in the 2D and 
3D systems. For example, the looped configuration in the 
2D system cannot undergo an out-of-plane buckling that 
would eliminate loops, but a 3D polymer will exhibit this 
behavior. Furthermore, in the 2D case it is much easier 
to form large defects which would be responsible for a 
larger decrease in extension. Therefore, the difference 
in the thermal-first or disorder-first ordering may be re- 
duced in a 3D system. But finally let us point out that 
the studies of the conformations of biopolymers are often 
performed in a 2D environment (e.g., see @, so our 
main findings should be instructive. 



Appendix 1: Proof of Eq. ([5]) 

Using the standard path integral methods [30| . for ar- 
bitrary function .F[{(^(s)}], we can write 

I V[<t>(s)]F[{<t>(s)}] cx lirn^ H I #^[{^}],(50) 

r r i 

/ V[<j>(s)]F[{<i>( S )}] (x lim J] / -^.F[{^}](,51) 
/ N—>oo / e 

where e = L/N, <fij = <fr[(j — l)e] is the discretized <fi(s), 
= <j)j — 4>j—\, and <f> in B and £ must be replaced by 
(<t>j+i — <t ) j)/ e - The Jacobian determinant, J = \d£/d<fi\, 
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of £'s with respect to </>'s is a constant, therefore, 



Appendix 2: A direct derivation of Eq. (|16[) 



JV-1 

2?[0(a)]^[{0( S )}] ex Jim -jp-j / #^[{^}[52) 

TV— »oo 6 / 
3=1 



Now taking J 7 = 73[{</>(s)}]e £ or e £ , from Eq. [[52]) we 
obtain Eq. 



/P[^( a )]B[{0( S )}]e- £ _ /P[^( S )]i?[{0( S )}]e- £ 



/23W*)]e- 



/P[0( S )]< 



(53) 



Intuitively, both sides in the above equation are averages 
over all possible configurations so they are expected to 
be equivalent. 



In this app endix, we will use the standard path integral 
method [3fJ to derive Eq. I]16p since it is useful. To 
account the more general case of sq ^ 0, we rewrite Tig 
as 



Ho = \ f *#( S ) - c(s)] 2 ds. 



(54) 



For large N, the path integral can be approximated as 



jv-i . 

z h ~ c n / d ^ ex p 

3=1 



N-l 
3=0 



/" 1 ™~ X 



where C is a constant, e = (L — sq)/N, <fij = cj)[so + (j — l)e], Kj — k[sq + (j — l)e] and c, = c[sq + (j — l)e] are 
discretized <f>(s), k(s) and c(s), respectively, and <^> in Tig is replaced by {4>j+\ — 4>j)/ e - Now using the identity 



( £ 2 , e -o(a;-zi) 2 -k(:E2-:E) 2 



a + b 



exp 



1/a+l/b 



xi - x 2 ) 



(56) 



we obtain 



c 



[ dfodfa ■ ■ ■ d0Ar_iexp \ V - 0j - Cj-e] 2 I 



•exp 



1 



2e/K + 2e/«i 



C 



•exp 



2?re / 2?re 



" 1 I 

2t/K l 1 + 2e/n 2 1 
1 



*o - (cq + ci)e] 



t>3«(p4 • • • d0Ar_iexp 



"3+i - W - c i 



h-<fo- (cq + ci + c 2 )e] 



= C exp < - 



2 i; 7 =o £ / k 3 



iV-1 



-*>-E 



3=0 



(57) 



where 1/k'i = l//to + and C is a new constant. Now let N — > oo, we obtain 



JE^ = C exp ■ 



2i4(L, s ) 



</>(L) - </>(s ) - / dsc(s) 

so 



(58) 



Normalizing the above equation we obtain C = 1/ y / 2irA(L, so), and recover Eq. 
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Appendix 3: Derivation of Eqs. (fl8)) -(f22 j) 

When force is small, expanding c~ £ about / = 0, we obtain 



P[0]e- £ «Z°(1 + /Q), 

B = ^-Jv^}B[m}e F J V{c}±-W({c}y 
— Bi — £?2, 

where 



7 7 l> 



V[c]W({c})(l-fQ)e- £o 



Q 



ds (cos <fi(s))° k , F=f coscj)ds, 



V^B^e 1 



V[c]W{{c})e 



-So 



w / V[<P\B[{<j>}]e 



-H 



» (B) K + f / <k'(£[{0( a )}] costs')),,, 
Jo 

B * = ^4o / mB[W)e F R^ I V^B^n, 

U = J V[c\W{{c})Qc- £o , 
and (••■)° denotes the ensemble average with energy £q 



(...)» = 4 / 2?[^(«)](...)e 



•70 
fe 



(59) 

(60) 
(61) 

(62) 
(63) 
(64) 

(65) 



Furthermore, 



n 



ds' 



ds 



G 

^Z 



H ° / 



g^e"^ / ds' 



V[c]W({c})e- £a[Ws)}] -^ [ 2?[0']cos[0'(s')]e~ £o[W ' (s)}1 
V{4>'] cos[^(s')] / P[c]T^({c})e~ £o[W(s)}I ^ o[W ' (s)}1 
cos[0'( s ')]e~ 1/2 /(f ds K^ 8 )"^ 8 )) 2 ] 
Jvm co S [</)'( s ')]e" 1/2/oids [ fc '( ' (s) "^ s) ) 2 ] 



7 70 



H 



t e - n ° / ds 



/D[^]e -1 / 2 -/c W>'(*)-7«] 2 
"W)/ 2 costs') +&)] 



(66) 



where 

fc(q + fc) 
fc ( s ) = — Zol - ' 7 ( s ) = 



7 («) = / 7(*)d«, G= v^c- 1 / 2 ^^ 2 ^ 2 , A'(s,s )= / d»/A/(«), (67) 



and we have used an expression similar to Eq. ([To]) in the last two lines in Eq. (f6"6"|) to transform the path integral 
into simple integral. j(s) is in general dependent on <f>(s) and it makes the expression complex. Now the Eq. (|63j) 
can be reduced into 



B 2 « X / ©MB[{0(*)}]e- w " 

TL 



ds l e -A\ s ',Q)/2 cos [ 7 ( s ') + ^ 

= // d S 'e-^'( s ''°)/ 2 ( J B[{0( S )}]cos[ 7 ( s ') + 0o]) K . 
Jo 

where </> = 0(0). Eqs. (f6"2")l and l[6"gj) are exactly the same as Eqs. (fT§)) and (|2"0"j) . 



(68) 
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Appendix 4: Calculations of the extension 

In this case B[{(f>(s)}} corresponds to cos0(s), and so from Eqs. (fT9|) - (|20p . we obtain 

X = (x)- 



= J I ds [ ds' 
>o Jo 



L pL 

ds (cos 4>{s)) — I ds (cos 4>(s)} t 



(cos(/>(s) cos^(s')) f 



C - S '' 2k ' (cos^s)] costs') + fa}) f 



(69) 



Note that Eq. (fT!)]) is valid only if s > Sq. Therefore, the integral for s' in Eq. ([74")) should be divided into two 
parts, one is from to s and the other is from s to L. When s > s', we have 



cos 4>(s) cos cj>(s')) K = / dcj)d(j)' P(4>, s; (f)' , s') cos 4> cos (f)' P((j)' , s';4> ,0) 



b' cos 2 ct>'P((t>', s; fa, 0) = - e - (s - s ')/ 2K [l + c- 2s ' /K cos(20 o ) 



and 



= e -(«-«')/2« 



cos 4>(s) costs') + fa]) K = I d4>d4>'P((j), s; <f>' , s') cos</>cos [k(</>' — fa)/a + <f> ] P(<f>' ', s'; fa, 0) 
-(s-s')/2 K 1 ,,,/«, )S ./ cos [«;(</>' -<^ )/q; + ^,] /'( <>'. s ':..„.() i 

cos(20 o ) + e 2s ' /c 



_J_ e -(s-s')/2K-(a+n) 2 6''/2KQ 



When s < s' , we obtain 

(cos0(s) cos <j>(s')) 

- p -(«'-«)/2« 



dcfydcf)' P((f>' , s'; 0, s) cos ^ cos (j>'P{<f>, s; fa, 0) 



cos 2 </)P(0, S ;</) O ,O) = -e-^- s )/ 2K 



1 + e- 2i '/ K cos(20 o ) , 

(cos0(s) costs') + ^o]) K = / dcj)d(j)'P((j)', s'; 4>, s) cos^cos [n(<f/ - fa)/a + fa] P(<f>, s; fa, 0) 

K(s'-,)/2a' I ,i, ,rn:<, ,COs[n{fa- fa)/a + fa]P(<>. s: ,, u . 0) 

cos(2</> ) + e 2s/c 



jj_ c -K(s'-s)/2Q 2 -(Q + K) 2 s/2fCQ 



It follows 



X = 2k/ 



L _ k(k + 3a) + k(k + a) c _ L/k _ Ana ^_ l/2k 
k + a k — a k — a 



k 2 cos(2fa)e~ 2L ^f 



X. 



3(2fc + a) (3k + a) 

X = 6k 2 (e L ' k - l) + ka (9e 2i / K - 16e 3L / 2K + 12e L / fe - 5) + a 2 (3e 2L / K - 8e 3i / 2K + 6e L / k - l) 
In the limit a — > k, X is still finite and 

X(a k) = 2k/ L-2k + (L + 2k)e- L ' k - -^fk 2 cos(20 o ) (l ~ 2e- L / fc + 2 C - 3i / fe - C - 4L / fe ) . 
On the other hand, in disorder- first system, 

X' = (%)'-{x)' f=0 = / ds {cos </>(s))' - / (coscf)(s)) K 



= f I ds / ds' (cos <j)(s)cos<j)(s' )} K — f 



ds (cos[4>(s)]} t 



2nf 



L + K -k (e- L ^ - 2)' + ^1 (e- 2i / K - & e - L ' K + 8e~ L ^ K - 3) cos(20 o ). 



When a — > 00, it is clearly that X — * X' , 



(70) 



(71) 



(72) 



(73) 
(74) 



(75) 



(76) 
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Appendix 5: Calculations of the orientational correlation function and the end-to-end Distance 

In this case, B[{<j)}} = t(s) • t(s') = cos[0(s) - (f>(s r )]. From Eqs. (fl8] ) -(|22 j) . we obtain 
(t(a)-t(a')) « (cos[0( S )-0( S ')]) K 

+/ f ds" [<«*#(«) - <f>(s')} cos0( S ")) K - e- s "/ 2k ' (cos[<f>(s) - </>(s')] cos[ 7 (s") + O ]> J • (77) 
The first term in above equation is simple, as 

/>OC 

(cos[0(a) - = / d^'P(0, S ;0 , , S , )cos(0-0 , )^'^ , ;0o,O) = e-( s - s ')/ 2K . (78) 



Again, due to that Eq. (JT6J) is valid only if s > ao, the integral for s" in Eq. (j77|) should be divided into several parts. 
If s > s' > a", we have 

(cos[0( S )-0( S ')]cos[ 7 ( S ") + ^o]) K 

d(pd<j)'d<j)"P{(j), s; s')P(</>', s'; 0", a") cos(> - <j/) cos[/t(</>" - 4>o)/a + <fo]P{<i>'\ a"; <fo, 0) 

/•OO 

= e - (s - s ' )/2K / #"cos[ K (0"-0o)/a + (/. o ]P(0",s";0o,O) =e- (5 '- s ' )/2K - ,i5 '" /2Q2 cos0o- (79) 



If s > s" > s', we obtain 

(cos[0( s )-0(s')]cos[ 7 (a ,, ) + <M) K 

d<t>d<t>' d4" P{4>, s; 0", s")P((f>", a"; <j)' , a') cos(<?(> - 0') cos[k(</>" - o )/a + (j> ]P((f>' , a'; (p , 0) 



= e-( s - s )/ 2k / d(j>'d(j>"P ((/)", a"; <p', s') cos(<£' - 0") cos[k(0" - <fo)/a + <fo] W, s'; <j) , 0) 

J — OO 

- L c -(s-s")/2n-( a+K ) 2 (s"-s')/2na 2 f l+c 2(s"-s')/ a \ / d(j)' COs[k((/>' - <j) ) / a + (j)o}P(<j)' , s' ; fa, 0) 

2 ^ ' J-oo 

= }_ c -s/2K+(3k+ a )s'/2ka-(3k+2a)s"/[2a(k+ a )] ^1 + e 2(*"-s')/«^ COS <£ . (80) 

If s" > a > a', we find 

(ca#(«) - cos[ 7 (a") + 0o]) K 

d<j)d<j)'d<j)"P((t)\ s"; <j>, a)P(0, s; a') cos(</> - (j)') cos[k(0" - o )/a + <fo]P{<l>' , s; <f>o, 0) 

/CO 
d<pd<p'dP(<j), s; <t>', a') cos(<?(> - 0') cos[k(0 - <^ )/a + 0o]P(</>', s'; cfo, 0) 
-OO 

1 / \ f°° 

= c -^"s)/2 a 2 -( a +^) 2 (s-s')/2^ f 1 + e 2(s-s')/ a \ / d4>'cos[K{(j)' -<t> )/a + (t) ]P{ct>',s';4> ,Q) 

* ^ ' J — oo 

= }_ e -(3k+ a )(s-s')/2k a - K s"/2a? ^ + e 2(W)/c^ COS0Q _ (81) 

Similarly, if s > s' > s", then 

(cos[0(«) - </»(*')] cos (<t>(s"))) K = e -( s - s ')/2«- s "/2K COS( ^ . (82 ) 

If a > a" > a', then 

(cos[0(s) - 0(a')] cos (#«"))>« = I e -( s - 3s '+ 3s ")/2« (i + c 2(s"-s')/^ CQS(j)o _ (83) 

If a" > a > a', then 

(cos[0(a) - 0(«')] cos (#«"))) „ - I e -(3 S -3,'+ s ")/2K ^ + e 2(.-«')/^ cos ^_ ( 84) 
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Combining Eqs. (177)) - ([84)) , we finally obtain for s > s' 

-S, 



(t( s ).t( s ')> « e-(*-°'^+ ' i ' a,S "' ,/ 



3(3fc + a)(k - a) 

S = -6fc(fc - a ) e -^+^)s+( a +3k)s']/2k a + ga ( 2k + a ) e -^/2K _ Qk ( 3k + a ) e -[2a S +(fe-a) S ']/2fea 



(85) 



+ (k - a)(3k + a) 



3 e -[(a+fc)L + (a+3fc)(s-s')]/2fca _ g e -[L+3(a-a')]/2« , 5g-s72« _ 3 e -(L+s'-s)/2K 



_|_2c _ ( 4 ' s_3 ' s ')/ 2k + 3 e -[(a+fc)L+(a-fe)(s-s')]/2fea 

(t(s) • t(s')) is also finite when a — > fc because in this case, 
(t(a)-t(a')) « e~ (s - s ' )/2K + KCOS(/.o/ 



6 2 



+e -[L+2( s - s ')]/fe + _ e (-4a+3a')/fc + e -[i+3(a-«')]/* + f -e"^ 

3 fc 



(86) 



In contrast, 



(t(s) ■ t(s'))' « (cos[0( s ) - 0(s')]) K + / / ds" (cos[0(s) - <f>(s')} cos^( s ")) K 

Jo 

-/<cos[0( S )-0( S ')]> K [ d s (cos[^s)]) K 
Jo 



-( s -a')/2« + y-y ^ cos [0( s ) _ 0( s ')] cos^(s")) K - 2 K /cos0 o e" ( ^ s ' )/2K (l - e~ L / 2K ) 

-(a-a')/2« + I K / COS O /'6e~ s ' /2K - 8e- s/2K + Q c -(L+s-s')/2n + 2c (3s'-±s)/2k 

3 V 

3e -[i+3(s-«')]/2K _ 3 e -(L+s'- s )/2 K ^j _ 



When ol — > oo, we obtain (t(s) • t(s')) -> (t(s) • t(s'))'. 
The end-to-end distance can be found by 

i? 2 = / ,/> •/ J.s'{t(>0 ■ tV)) 2/ <Is I ds' {cos[(j>(s) - 0(s')]) 

4k 2 cos 0o / 



= 4kL 



2k , 
l- T (l-«. 



(88) 



9(fc - a) 2 (fc + a)(2fc + a) (3k + a) 3 ' 
y = 18(k-a) 2 (k + a)(2k + a)(3k + a) 2 L~9k(k-a) 2 (3k + a) 2 (4:k 2 + lla 2 + 22ka) 

-16a(2fc + a) 2 [3(fc - a)(k + a) (3k + a)L + 4ka(5k + a) (3k - a)]e" i/2K + 18ka{k + af(k - a ^ e -(i/k+2/ a )L 
+18k(k + a) 3 (2k + a)(3A + a) 2 e - L / fc - ka(k - a) 2 (2k + a)(3k + a) 2 e" 2L/K . (89) 



when a — > k. 



R 2 = 2kL 



(l-e-*) 



l--[l-e- 



+ 



k cos ipo f 
18 



36kL - lllfc 2 + (18i 2 + 78fcL + 109fc 2 )e~ L/fe + 3fc 3 e" 3i/fc - k 2 e~ 4L/k 

(90) 



On the other hand 



R' 2 = AkL 



2k ( 

1 ~ T ( 37 



-k 2 cos ^o/ 18L - 99k + 16(3L + 4k)c- l / 2k + 36kc- l / k - kc- 2L / k . (91) 



When a — > oo, we obtain R 2 — > i?'' 
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